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,_i ■ The complete analysis of a model with three quartic coupling constants associated 

with an 0(2iV)-symmetric, a cubic, and a tetragonal interactions is carried out within 
the three-loop approximation of the renormalizat ion-group (RG) approach in D = 
4 — 2e dimensions. Perturbation expansions for RG functions are calculated using 
dimensional regularization and the minimal subtraction (MS) scheme. It is shown 
that for N > 2 the model does possess a stable fixed point in three dimensional 
space of coupling constants, in accordance with predictions made earlier on the base 
of the lower-order approximations. Numerical estimate for critical (marginal) value of 
the order parameter dimensionality N c is given using Pade-Borel summation of the 
corresponding e-expansion series obtained. It is observed that two-fold degeneracy of 
the eigenvalue exponents in the one-loop approximation for the unique stable fixed 
point leads to the substantial decrease of the accuracy expected within three loops 
and may cause powers of \/e to appear in the expansions. The critical exponents 7 
and 77 are calculated for all fixed points up to e 3 and e 4 , respectively, and processed 
by the Borel summation method modified with a conformal mapping. For the unique 
stable fixed point the magnetic susceptibility exponent 7 for N = 2 is found to differ 
in third order in e from that of an 0(4)-symmetric point. Qualitative comparison of 
the results given by e-expansion, three-dimensional RG analysis, non-perturbative RG 
arguments, and experimental data is performed. 

PACS numbers: 64.60.Ak, 64.60.Fr, 75.40. Cx, 75.50.Ee 
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1 Introduction 

There are numerous complicated models with more than two independent quartic cou- 
pling constants. They describe phase transitions in a variety of systems and are actively 
studied within the e-expansion as well as by the field-theoretical renormalization group 
method in three dimensions [1-10]. Critical fluctuations in anisotropic systems with 
several quartic coupling constants are known to destroy, as usual, continuous transi- 
tions converting them into first-order ones. This fact, however, has not been strictly 
proved. On the contrary, it is possible to construct models with a large number of 



coupling constants whose RG equations have stable fixed points |lT[. It means that 
the presence of three and more coupling constants in the Landau-Wilson Hamiltonian 
does not forbid continuous phase transitions in the system. Nevertheless, complicated 
models with stable fixed points are quite rare. One of such models describing certain 
antiferromagnetic phase transitions and the structural transition in Nb02 crystal will 
be studied in the paper. 

We consider the critical behavior of a model given by the fluctuation Landau- Wilson 
Hamiltonian with three quartic interaction terms: 






H = d D x 

1 

4! 



Here (pf, i = 1, . . . , N, a = 1,2, is the real 2A r -component order parameter field in 
D = 4 — 2e dimensions and m , u , v , z are the bare mass and coupling constants, 
respectively. The squared bare mass m^ can be thought of as proportional to the devi- 
ation from the mean-field transition point (line). The field (p® is regarded as consisting 
of two sets of components, even and odd, each of them may be considered as a real 
iV-component vector. Tensors G±, G2, and G3 in the Hamiltonian (^ corresponding 
to isotropic, cubic, and tetragonal interactions have the following symmetrized form: 

Gi gfcf = ~ (^""Mw + S a ^SikSji + S au ^Suhj) , 

G 2 fkT = 5 af3 5 a »5 a »5 tl 5 tk 5 th (2) 

G 3 gT = - (5 la 5^5 2fl 5 2u + 5 la 5^5 2(S 5 2u + 5 la 5 lu 5 2 » 5 2f3 

+5 2a 5 w 5^5 lu + 5 2a 5 2 »5 1/3 5 lu + S 2a S 2u S ll "S 10 )6 ij S ik S u . 
When u is equal to zero, the Hamiltonian (Til) describes iV identical non-interacting 



anisotropic XY models [12 , while for z = it turns into the Hamiltonian of the 



well-known hypercubic model. 



The Hamiltonian ([D) governs the critical thermodynamics in a number of interesting 
physical systems. So, for example, when N = 2 it describes the structural phase tran- 
sition in Nb0 2 crystal and, for v — z, the antiferromagnetic transitions in TbAu2 and 
DyC 2 - Another physically important case iV = 3 is relevant to the antiferromagnetic 
phase transition in K 2 ItCIq crystal and, for v — z, to those in TbD 2 and Nd [[!} |3|. 
The detailed analysis of these systems along the line of the Landau phenomenological 
theory can be found in Jl|, |2|, [13|] with references to the experimental works therein. 

For the first time the renormalization group analysis of the model ([I]) was performed 
to second order in e by Mukamel and Krinsky in Refs. [|l], §, |!|. On this ground, it 
was shown that the 2iV-component real anisotropic model (|I|) has a unique (three- 
dimensionally) stable fixed point for each N > 2. The corresponding critical exponents 
were recorded and for n = 2N = 4 they were found to coincide with those of the 
Heisenberg fixed point. On the other hand, the critical behavior of this model was 
studied within the two-loop approximation by the alternative RG approach in three 
dimensions |14]]. The calculations made provided the same qualitative predictions, 
although for the physically interesting cases N = 2 and N = 3 the critical exponents 
were found to be numerically close to those of the 3D XY model rather than the 
Heisenberg ones. 

It is well known, however, that low-order approximations lead to rather crude 
quantitative and, sometimes, contradictory qualitative results, especially for systems 
with nontrivial symmetry (see, for instance, Refs. [|7], |^, |T3[). To make more definite 
conclusion concerning the unique fixed point stability and obtain more accurate values 
of the critical exponents one has to consider long enough perturbation theory series. 
Such series are known to have the zero radius of convergence and therefore are, at 
best, asymptotic. To extract reliable information from them a proper resummation 
procedure must be applied. Recently such work for the model under consideration was 
done within the field-theoretical RG approach in three dimensions, where the three- 
loop expansions for /^-functions and critical exponents were calculated for arbitrary iV 
[ |16|j . Using the generalized Pade-Borel resummation technique, the coordinates of all 
the fixed points were found. It was shown that the unique stable fixed point did exist 
on the three-dimensional RG flow diagram when N >2. 

It should be noted that, assuming v — z the model ([]]) formally turns into that 
with generalized cubic anisotropy and the complex order parameter field. The latter is 
a specific case (m = 2) of the well-known rrm-component model. The critical behavior 
of this model was investigated in Refs. |T7|, |IBfl . Two- and three-loop calculations done 
for the case m = 2, n > 2 predict stability of the mixed fixed point, the analog of the 
unique stable fixed point of the model ([I]). 

At the same time, there are general non-perturbative arguments in favor of the 
unique stable fixed point should not be in the physical space although its existence 



is not forbidden at D > 3 ||19| . According to those considerations the only three- 
dimensionally stable fixed point may be the Bose one and it is that point which governs 
the critical thermodynamics in the phase transitions mentioned. The point is that the 
model ([!]) describes N interacting Bose systems when v — z. As was shown by J. Sak 
[p0|| , the interaction term can be represented as the product of the energy operators 
of various two-component subsystems. It was also found that one of the eigenvalue 
exponents characterizing the evolution of this term under the renormalization group in 
a neighborhood of the Bose fixed point is proportional to the specific heat exponent a. 
Since a is believed to be negative at this point (that is confirmed by highly precise up- 
to-date mesurements of the specific heat exponent of liquid Helium |?TJ including those 



in outer space [22] and the high-loop RG computations carried out for the simple 0(n)- 



symmetric model in three dimensions E31, E3J) the interaction is irrelevant. Therefore 



the Bose fixed point should be stable in three dimensions. 

Renormalization group approach, however, when directly applied to the model (|l|) 
and to the relative win-component one, has not still confirmed that non-perturbative 
conclusion. On the contrary, all calculations performed up to now indicate existence 
of the unique stable fixed point in the physical space, while the Bose point appears to 
be three-dimensionally unstable [1-3,14,16-18]. This may be a consequence of rather 
crude approximations used, and the higher order being taken into account the closer 
the perturbative results could be to the precise ones. So, the aim of the paper is to 
investigate the critical behavior of the three coupling constants model (ffl) in the next, 
third order in e and verify compatibility of predictions given by the e-expansion method 
with predictions based on the other techniques. 

The main results of our study to be discussed below are as follows. 

• The /3-functions of the record length for the model ([!]) are obtained by the e- 
expansion method. To calculate tensor convolutions associated with the Feyn- 
man's graphs an algorithm was developed and a specially designed computer 
application package was written. 

• Coordinates of all fixed points and their eigenvalue exponents are calculated in 
general form within the three-loop approximation. The problem of stability of 
the fixed points is analyzed. The unique fixed point rather than the Bose one is 
found to be three-dimensionally stable in the frame of given approximation. Nu- 
merical estimate of the critical dimensionality of the order parameter N c obtained 
confirms this conclusion. 

• It is observed that one-loop degeneracy of the eigenvalue exponents of the unique 
stable fixed point leads to a certain complication in calculating their e-expansion 
series. This problem is investigated in detail. It is shown that such a degeneracy 



substantially reduces the accuracy expected from given approximation and may 
result in appearance of the powers of y/e in corresponding series. 

• Perturbation series for the critical exponents 7 and r\ are expanded to e 3 and e 4 , 
respectively. For N = 2 the magnetic susceptibility exponent series of the unique 
stable fixed point and the 0(4)-symmetric point are found to be different (up to 
second order in e they exactly coincide ||). The numerical values of the critical 
exponents are estimated by resumming the series using the Borel transformation 
with a conformal mapping. 

The results of our investigation are discussed in comparison with conclusions given by 
other theoretical approaches and experimental data. 

The set up of the paper is as follows. In Section [| the renormalization scheme is 
formulated and three-loop expansions for the /^-functions and critical exponents are 
presented. Specific symmetry properties of the initial Hamiltonian (|l]) are revealed and 
used as criteria of correctness of the equations deduced. In Section || RG expansions 
for coordinates of the fixed points are written out for arbitrary N and the problem of 
their stability is studied. The numerical estimate of the critical dimensionality N c , at 
which the topology of flow diagrams changes, is obtained in Section [| by resummation 
of its series using Pade-Borel method. The RG expansions of the critical exponents 
for the physically interesting cases N = 2, N = 3 and their numerical estimates are 
given therein. Conclusion is devoted to discussions of the results of the investigation. 
The paper has two appendices. Appendix A contains e-expansions for the eigenvalue 
exponents of the fixed points for arbitrary N. In Appendix B we analyze the problem 
of degeneracy of the eigenvalue exponents of the unique stable fixed point and its 
implications. 



2 RG expansions and symmetries 

To calculate the /3-functions and critical exponents normalizing conditions must be 

(2) 
imposed on renormalized one-particle irreducible inverse Green's functions T R and 

vertices T^ R given by corresponding Feynman's diagrams. Within the massless theory 



they are normalized in a conventional way: 

Fr\{p};h,u,v,z) 

ru R ({phn,u,v,z) 

Fzr({p}]H,u,v,z) 
with one more condition on the if 2 insertion 

r ( R 2 \{p},{q};H,u,v,z) 




1 

fi e z 



(3) 
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(4) 



Here m, u, v, and z are the renormalized mass and dimensionless coupling constants, 
with fj, being an arbitrary mass parameter introduced for dimensional regularization. 

The vertices r^ 4 -*, ly-*, T^> are connected with the vertex function without external 
lines normalized in the following way: 



-,(4) afinv 
ijkl 



r (4) ft a/3fMi> , p(4) (~t aj3au , r (4) (~t 
L u ' ^lijkl "r L v ' (j 2ijkl ' L z ' *** 



r 3 ijkl ' 



From renormalization conditions @ and (|J) the expansions for the renormalization 
constants Z^, Z w Z v , Z z , and Z^p. may be obtained. These constants relate the bare 
mass txiq and three coupling constants Uq, Vq, z of the Hamiltonian ([!]) to the corre- 
sponding physical parameters: 



m r 



Z<p2 2 

—=—m 



Z m m 



Z l 



Z l 



Z l 



(5) 



Thus, with relations (|5|) taken into account, the /3-functions and critical exponents can 
be calculated via the formulas 



G>lnitp a _i_ d In up o _i_ d In up o 
du t Ju "I" dv rvT dz t Jz 

dlnvg a I d In vp a , dlnvp n 
g u > Ju ^ dv Vv -r g 2 ^z 

d In z /o _|_ dlnzo/o , d In z a 
du l Ju ^ dv l Jv ^~ dz ^ z 



-e, 

-e, 



(6) 



2^/3 u + 2^p v + 2^p z 

.dlnZo „ „31nX„2 „ „dln] 



(7) 



rj(u,v,z) = 

m(%v,z) = 2"-^f3 u + 2"-^{3 v + 2"-^f3 z 

-2— g = u,v, z. The critical exponents r\ and 772 are found by substituting 



where /? g = 9|lnM| 

zeros of the /^-functions into expressions (f 

well known scaling relation 7 -1 = 1 + ^-. 



). The critical exponent 7 is given by the 



The contribution of a Feynman's graph into an RG-function comprises three fac- 
tors: the combinatorial coefficient, the result of tensor convolution and the integral 
value associated to the diagram. The combinatorial factors, and the values of integrals 



are known from Ref. [25], while evaluating tensor convolutions for vertex and mass 
diagrams is the problem to be solved. To do it we have developed a computer appli- 
cation package written in PASCAL. The algorithm is based upon two quite natural 
assumptions: 

1. Tensor convolution algebra is closed, i.e. each monomial G^ * ... * Gi l+1 con- 
tributing to a vertex function is a linear combination of the basic tensors G±, G2, 
and G3: 

G n *...*G ll+1 = a(N)G 1 + b(N)G 2 + c(N)G 3 . (8) 

2. Dependence of the coefficients a(N), b(N), and c(N) upon N is of polynomial 
character. The degree of the polynomials does not exceed the number / of loops 
in the Feynman's graph. 

The first condition means that one has no new interactions generated in the model 
([I]). The second proposition becomes evident upon analyzing the particular form of 
the tensors Gj. 

Since a polynomial of degree I is fully determined by its values in / + 1 different 
points, it is sufficient to compute convolutions consecutively assuming N = 2, . . .1 + 2 
(the reason to start from 2 is linear dependence between G{ when N = 1). In order 
to evaluate three indeterminates a(N), b(N), and c(N) we compare both sides of 
expression (f|), having assigned values (iili); (1111)) an d (1111) to the multi-index Qt^ ). 
It provides a non-degenerate system of linear equations whose 3 x 3-matrix does not 
depend on N. From this system the coefficients of decomposition (|8|) are found. Similar 
procedure was applied to the mass diagrams. The results of our computations recover 
those achieved within the four-loop approximation for simple 0(n)-symmetric model 

After some work, we obtain the expressions for the RG-functions within the three- 
loop approximation (Fisher's exponent r\ is calculated up to four loops) using dimen- 
sional regularization |26| and the MS scheme |27j : 



(3 U = eu-u 2 - 2{I } + 4) ( Quv + 2uz ) + 

4(JV ] f4) . [l2w 3 (3iV + 7) + I32u 2 v + Uu 2 z + 30m; 2 + 10uz 2 
16(jv 1 +4)3 [4m 4 (48C(3)(5AT + 11) + 33iV 2 + AQ1N + 740) + 
12u 3 w(384C(3) + 79iV + 659) + 4u 3 z(384C(3) + 79iV + 659) + 
18uV(96C(3) +N + 321) + 1380u 2 uz + 2m 2 z 2 (288C(3) + 3iV+ 
733) + 1512m; 3 + 18uv 2 z + 504uvz 2 + 222uz^ 



Pv 



ev — 



(12uv + 9v 2 + z 2 ) + 



4(AT+4) 2 
8^ 



rp [4u 2 v{5N + 41) + 276uv 2 + 20uvz + 24uz 2 + 102w 3 + 10vz 2 + 



16(^+4)3 



8M 3 w(96C(3)(iV + 7) - 13iV 2 + 184iV + 821) + 



l8wV(768C(3) + l7iV + 975) + 12m 2 ^(96C(3) - 137V + 154) + 
2uV(576C(3) + 43iV + 667) + l08m; 3 (96C(3) + I3l) + 30Quv 2 z+ 
12m^ 2 (96C(3) + 187) + 2wz 3 (384C(3) + 395) + 27v 4 (96C(3) + 145) + 
I62v 2 z 2 + 8^ 3 (48C(3) + lOl) + 3^ 4 (32C(3) + 17) 1 



(9) 



(3 Z 



7 



-i 



ez 



2(iV + 4) 



'12uz + Qvz + Az 2 ) + 



4(iV + 4) 2 
1 



4w 2 z(5iV + 41) + 204uvz + 116uz 2 + 30v 2 z + 72vz 2 + 18z s 



16(iV + 4) 3 



8w 3 z(96C(3)(iV + 7) - 13iV 2 + 184iV + 821) + 

12w 2 ^(864C(3) + AN + 1129) + 4uV(1440C(3) + 47iV + 1796) + 
18w 2 z(192C(3) + 391) + 72w^ 2 (96C(3) + 103) + 2m/(960C(3) + 1517) 
1512v 3 z + 36vV(48C(3) + 35) + 72v* 3 (16C(3) + 25) + 4/(48C(3) + 91) 



1- 



1 



2(iV + 4) 
1 



(2u(N 



Sv + z) + 



2(iV + 4) 2 
1 



6u 2 (N + 1) + 18uv + 6uz + 9v 2 + 3z z 



\2u 3 {N + 1)(11JV + 39) + hlu 2 v{llN + 39) + 



16(iV + 4) 3 

\Su 2 z(llN + 39) + 6w 2 (5iV + 398) + hQAuvz + 

2wz 2 (5iV + 304) + 80k; 3 + 15t» 2 ^ + 267uz 2 + 117/ 



(10) 



2{N 



-—-{2u 2 {N + 1) + Guv + 2uz + 3v 2 + z 2 ) 



4u 3 (N + 1)(N + 4) + 18u 2 v{N + 4) + 6u 2 z(N + 4) + 



3(iV + 4) 3 

Muv 2 + l%uvz + 2luz 2 + 27w 3 + 9vz 2 + 4z z 



:ir 



40u 4 (iV + l)(iV 2 -9N- 25) + 240u%(N 2 - 9N - 25) -} 



32(iV + 4) 4 

80w 3 ^(iV 2 - 9JV - 25) - 180uV(iV + 58) + 360u 2 vz(N 
180v? z 2 (N + 14) - 7020mw 3 - 180uv 2 z - 2340uvz 2 - 
1020uz 3 - 1755v 4 - 90v 2 z 2 - 720vz 3 - 75/ 



are in 



where ( is the Riemann (^-function: ((3) = 1.20206. Expressions @ - (ljj 
accordance with those obtained earlier in Ref. ||, where corresponding calculations 
for RG functions were carried out to e 2 , and, assuming v =- z = and iV = §, with 



results of Ref. ||28| , in which the critical exponents of the well-known 0(n)-symmetric 
model were found up to e 4 . If u = and t> = 2, the right-hand side of the second 
(third) equation @ goes over into the /^-function for Bose-like systems, the coupling 
constants being normalized properly. The latter, obviously, coincides with that of the 
0(n)-symmetric model when n = 2N = 2. 

In conclusion of this section, let us formulate a criterion to check the correctness of 
the results obtained. It relies on a specific symmetry property of the Hamiltonian (|l|) 
of the system under consideration [Q. It occurs that transformation 

<f2N-l — > 4| (f2N-l + <P2N), 

(12) 

¥>2N —> TyJ (<P2N-1 — <P2n), 

combined with substitution of quartic coupling constants 

1 1 

u^u, v->~(v + z), z-*~(3v-z) (13) 

does not change the structure of the Hamiltonian itself. 

Similar situation takes place for N — 1 and z = in the case of another field 
transformation 

(14) 

<^2 -> ^(<fl~<f2), 

which does not affect the Hamiltonian resulting only in the following replacement of u 

and v : 

3 

u — > u H — v, v — > — f . (15) 

It is well known that the RG equations should be invariant with respect to any 
transformation conserving the structure of the Hamiltonian [25|. It means that for 



every N, in the case of symmetry (|i3|), functions (3 U , (3 V , and (3 Z should obey special 
relations which may be readily written down: 

(3 u (u, v, z) = (3 u \u, \ (v + z), \ (3v-z) 
f3 v (u, v, z)+(3 z (u, v, z) = 2(3 v (u, ~ (v + z), \ (3v - z) ) , (16) 

3f3 v (u, v, z) - (3 z (u, v, z) = 2f3 z (u, ~ (v + z), \ (3v - z) 

For N = 1 and z = the other symmetry (fL5|) results in 

/3 u (u, v, 0) + ~f3 v (u, v, 0) = /3 u (u+~v, -v,0 

P v (u, v, 0) = -(3 v iu+\v, — v, 

At last, the critical exponents are invariant under the transformations dl3|) and (jH 



(17) 



So, the first symmetry gives 

7~ 1 (w, v, z) = 7 _1 fM, \ (v + z), \ (3v-z) 
r](u, v, z) = r)(u, \ (v + z), \ (3v — z) 

Similar relations should take place in the case of the symmetry dl5|) . It can be easily 
verified that conditions (|X6|) — (p~8|) are satisfied indeed. 

Symmetries of the initial Hamiltonian like those described above (such symmetries 
do not exist always and to find them requires certain efforts) play, in some cases, an ex- 
tremely important role. Namely, the series being obtained within the field-theoretical 
RG approach in 3D are necessarily processed with the use of some resummation pro- 
cedure (e.g. Pade, Pade-Borel, Chisholm-Borel etc.), and satisfaction of the numerical 
results to the exact symmetry relations serves as a criterion to estimate the accuracy 
expected from the approximation scheme employed f7[. 



3 Fixed points and stability 

Two critical exponents 7 and 77 are known to completely specify the critical behavior of 
a system [3(J. They are determined from RG functions by going to the infrared-stable 



fixed points g c = (u c , v c , z c ), which are found as zeros of the /3-functions in the form of 
series in powers of e: 



9c = 9c(e) = J2dke k . 
fc=i 



There exist eight fixed points in the model under consideration []3], |nj, one of them 
(Gaussian) is trivial: 

1. Gaussian fixed point 

u c = v c = z c = 0. 

2. 0(2N)— symmetric or Heisenberg fixed point 

_ , 3(3AT+7) 2 / 12C(3)(5AT+11) , 33iV 3 -55JV 2 -440JV-568 ^ r 3 

" c ~~ fc "r (AT+4)2 fc V (^+4)3 " 1 " 4(Af+4)4 J fc > 

V c = Zc = 0. 

3. Ising fixed point 



U c = Z c = 0, 

2(iV+4) e _,_ 68(iV+4) c 2 j_ / 709 (N +4) _ 32(JV+4) , 



.. . , .) , 68(iV+4) 2 I / 709(7V+4) 32(JV+4) >/ x\ 3 
" CH 2I3 - C + V 6561 81 U^JJ 6 • 



Cubic fixed point 

u c = ^e + |±|(l-2iV)(19iV-53) e 2 + 

( 4 g±£C(3)(8iV 3 - 12iV 2 -7N + 7)- 

^g^ (391(W 4 + 41971iV 3 - 114987iV 2 + 9016(W - 22472))e 3 , 
vc = ^^(iV-2) e + g±|)(2iV-l)(17iV 2 + 55iV-106)e 2 + 

(-^^C(3)(2iV 4 + 8iV 3 - 10iV 2 -9N + 7) + 

13 ^ 5 (1418iV 5 + 11713iV 4 + 90281iV 3 - 247414iV 2 + 187528iV- 



44944) je 3 , 
z c = 0. 
5. Bose fixed point 

u c = 0, 

v c = ^ eH -^v±4l e2H _^(io3-384C(3))e 



^ = A[±4 e+ 6(iv±4l e2 + ^ ( i03-384C(3))e 3 . 



5 ' 25 ' 1250 

+4 , 6(JV+4) a , JV+4 

5 ~ l ~ 25 ^ 1250 

VZ— cubic fixed point 

u c = 0, 

Wc = ^e+ 3 i^ e 2 + ^(709-2592C(3))e 3 , 



9 243 13122 

^ = ^e + ^%±^e 2 + f±f(709-2592C(3))e 3 . 
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7. I-tetragonal fixed point 

u c = ^e + !±t(l-2iV)(19iV-53)e 2 + 

(^#C(3)(8iV 3 - 12iV 2 -7N + 7)- 

gg^5 (3910AT 4 + 41971iV 3 - 114987iV 2 + 90160iV - 22472))e 3 , 
v c = ^(iV-2)e+^(2iV-l)(17iV 2 + 55iV-106)e 2 - 

(S#C(3)(2iV 4 + 8iV 3 _ 10iV 2 _ 9N + 7) _ _g±4_ (1418iV 5 + 

11713iV 4 + 90281iV 3 - 247414iV 2 + 187528iV - 44944))e 3 , 
z c = ^(iV-2)e+g±|(2iV-l)(17iV 2 + 55iV-106)e 2 - 

(^#C(3)(2iV 4 + 8iV 3 - 10iV 2 - W + 7) - g^(1418iV 5 + 
11713iV 4 + 90281iV 3 - 247414iV 2 + 187528iV - 44944))e 3 . 

8. II-tetragonal fixed point 

u c = (^ y 6+ Ii §| F (70iV 2 -205iV + 139)e 2 + 
(|gS 4 lC(3)(64iV 3 - 188iV 2 + 1517V - 23) + 
-4-5^)5 (6370iV 4 + 24149 A^ 3 - 144719 A^ 2 + 197208iV - 83256) 
jg&;{N - 2)6 + ^J^ (3(W 3 + 25iV 2 - 217iV + 166)e 2 - 



v, = 



4(4-5JV1 8 lu " u '' : ~ !J U7 " x ~~ ' ! u J '"'- * ' i-v<-uo.\ — ^."••)„-ni)ie 

JV+4_ 
4 

(f|^4yC(3)(8iV 4 + 16iV 3 - 88iV 2 + 75iV - 9) - i ^ F (1030iV 5 + 

e 3 , 



2751iV 4 + 46033iV 3 - 207590iV 2 + 267336iV - 109808) 
$±(N- 2)e + j^_(30N 3 + 25iV 2 - 217iV + 166)e 2 - 
(Sfc4yC(3)(8iV 4 + 16iV 3 - 88iV 2 + 75iV - 9) - i J^(1030iV 5 + 
2751iV 4 + 46033iV 3 - 207590iV 2 + 267336iV - 109808)^ 



From these expressions it is seen that for the physically interesting case N = 2 the 
coordinates of the fixed points 2, 4, 7, and 8 coincide in the one-loop approximation, 
i.e. the Heisenberg point u c = e, v c = z c = is four-fold degenerate. Such strong 
degeneracy is occasional, however, and lifted out in higher orders of the perturbation 
theory. So, the two-loop approximation splits those points apart. This situation is 
typical for a number of complicated models (see, for example, Refs. |7|, ||). 

One can also notice from the above list that the Heisenberg and the Ising fixed 
points coincide at iV = |, while for the components of the cubic and the Ising fixed 
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points the relation — v c c = v* holds at N — 1. With A" increasing, the cubic fixed point 
approaches the Heisenberg point from below and crosses it at N — N c , changing the 
sign of its f-coordinate. Further, when iV — > oo it moves towards the Ising point. 
Note, that the II-tetragonal fixed point is getting close to the Bose one when N grows. 
Such a behavior of the fixed points is in accordance with results obtained within the 



RG analysis in three dimensions [14, [1 



Since the symmetry transformations (|i3|) , ([15]) do not affect the form of the RG 



equations, they can only rearrange the fixed points. This observation may be used as 
an additional criterion for verification of our results. For example, points 1, 2, 5, and 
8 stay untouched under transformation (|IB|), while points 3 and 4 turn into 6 and 7, 
respectively (and vice versa). 

Now let us discuss the character of the stability of the fixed points found. It is 
known to be determined by the signs of the eigenvalues Ai, A2, and A3 of the matrix 



M l0 



d(5 u d(3 u dj3 u 



S0 Z d(J z d(3 z 
du dv dz 



evaluated at u = u c , v = v c , and z = z c . If the real parts of all the eigenvalue expo- 
nents are negative, the corresponding fixed point is infrared stable in three dimensional 
(u, v, z)-space. Besides, the "saddle-knot" type fixed points may occur on the phase 
diagram, provided their eigenvalue exponents are of opposite signs. General expres- 
sions for the eigenvalue exponents are written out for arbitrary N in Appendix A. For 
the interesting cases N = 2 and A^ = 3 relevant to the substances of concern they are 
presented in Table 1 and Table 2, where e = | corresponds to the physical case. 

It is seen from the tables that the Ising point has single negative eigenvalue, there- 
fore it is stable only on the f-axis. The Heisenberg point is stable on the axis too if 
A^ > N c , becoming stable within the plane (u, v) for A^ < N c . The cubic fixed point 
has the critical behavior opposite to that of the Heisenberg one; they interchange their 
stability at A^ = N c . The Bose point is stable within the plane (v,z), being of the 
"saddle-knot" type in the three-parameter space. Note that the eigenvalue exponents 
of points 3 and 6 are the same as well as those of 4 and 7. It is a consequence of the 
symmetry (|T^) of the initial Hamiltonian. 
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Table 1: Eigenvalue exponents for N = 2 to third order in e. 



No 


Type of fixed point 


Eigenvalues 


1 


Gaussian 


A u = X v = \ z = e 


2 


Heisenberg 


\ _ , , 13 r 2 84<(3)+65 3 

a u — t -r 12 t 3 g t 

\ - \ - l r 2 1 i-^r 3 
At, — A z — 3 e + 6 e 


3 


Ising 


X _ \ _ 1, 38 r 2 | 2592C(3)-937 3 

a u a 2 3 t gl t -r 21g7 

\ _ , , 34 r 2 2592<(3)+1603 3 

A v t -r 27 t 72g 


4 


Cubic 


X _ r | 13 r 2 84C(3)+65 3 

Ai - e + 12 e 36 e 

\ _ l r 2 | 15C(3)+1 3 

A 2 - 3 e + 18 e 

> _ l r 2 15C(3)-7 3 
A z — 3 t lg t 


5 


Bose 


\ _ 1, 14 r 2 | 768C(3)-311 3 
Au 5 fc 25 fc ^ 625 
X _ l r | 2 9 768C(3)+29 3 
Al — 5 C + 5 C 625 6 
X _ , I 6 r 2 384<(3)+257 3 
A 2 - e +5 C 125 6 


6 


^Z-cubic 


\ _ \ _ 1, 38,2 | 2592C(3)-937 3 

y\ u ^i 3 t gl t -r 21g7 

\ _ . i 34^2 2592C(3)+1603^ 3 

A 2 - e + 27 e 729 e 


7 


I-tetragonal 


\ _ 1,2 15C(3)-7 3 

Ai - 3 € 18 e 

\ _ l r 2 , 15C(3)+1 3 

A 2 - 3 e + lg e 

X _ , | 13 r 2 84C(3)+65 3 

A 3 - e + 12 e 30 e 


8 


II-tetragonal 


Ai = A 2 = — 3 e 

\ _ r | 13 r 2 84C(3)+65 3 

A 3 - e + 12 e 36 e 
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The most intriguing is the II-tetragonal fixed point proving to be absolutely stable 
in 3D, as it follows from the tables. Obviously, simple resummation procedures, such 
as Pade and Pade-Borel methods, applied to A's do not dismiss this conclusion. The 
presence of such a stable point is extremely important. It implies that the critical fluc- 
tuations do not destroy the second-order phase transitions, at least, if the anisotropy 
of the initial Hamiltonian is not too strong. Since the stable fixed point is located on 
the plane v = z it is certainly relevant to the critical behavior of TbAu2, D1/C2, TbD2, 
and Nd. 

Let us note, that the e-expansions of the eigenvalue exponents Ai and A 2 for the 
II-tetragonal point substantially differ from the others. Namely, their series prove to 
be shorter by one order (see Table 1 and Table 2). This phenomenon originates from 
the two-fold degeneracy of the roots of the characteristic polynomial in the one-loop 
approximation. As a consequence, the eigenvalue exponents should be expanded in y/e 
rather than e. It can be shown, however, that for almost all N non-integer powers 
drop from Ai and A2 for the eighth fixed point in every order of the perturbation 
theory. For the special case N = 2, significant from the physical viewpoint, one cannot 
make such a statement within three-loop approximation. To answer that question one 
should take into account at least four-loop contributions. Apart from whether or not 
non-integer powers of e appear in the expansions, one-loop degeneracy of A's results 
in reduction of information available from a given approximation. So, assuming that 
sfe will not appear in the series for N = 2, evaluating coefficients of A's in third order 
in e would require accounting five-loop terms. To understand the structure of the 
eigenvalue exponent series with one-loop degeneracy, we conduct detailed analysis of 
the problem in Appendix B. 



4 Marginal dimensionality and critical exponents 

We have shown in the previous section that for the physically interesting cases N = 2 
and N = 3 the II-tetragonal fixed point is three-dimensionally stable in 3D. The 
question may be put forward whether this point is stable for all N. To answer it the 
critical dimensionality of the order parameter N c needs to be calculated. It separates 
two different regimes of critical behavior of the model. When N > N c the II-tetragonal 
rather than the Bose fixed point is three-dimensionally stable in 3D. At N = N c they 
interchange their stability so that when N < N c the stable fixed point is the Bose one. 
The e-expansion for N c can be found from the condition v c = z c = imposed on the 
coordinates of the eighth fixed point (see Sec. |3|). Three-loop approximation gives 

N c = 2-2e + ^(6C(3)-l)e 2 + 0(e 3 ), (19) 
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Table 2: Eigenvalue exponents for N = 3 to third order in e. 



No 


Type of fixed point 


Eigenvalues 


1 


Gaussian 


A« = 


A^ — A z — e 


2 


Heisenberg 


A« — 

\ v = 


. 48 2 2(2184C(3)+1931) 3 
^ 49 2401 
\ __ 1 | 104 2 2(5208C(3)-2311) 3 
z 7 fc ' 343 16807 


3 


Ising 


A« — 
A v = 


\ __ 1, 38 2 | 2592C(3)-937 3 

A 2 gt gl t "I" 21g7 ! 

, 34 2 2592C(3)+1603 3 
"^ 27 729 


4 


Cubic 


Ai = 

A 2 = 


1 250 2 246888C (3) +165287 3 

t ' r 243 118098 
1 250 2 , 5(36936C(3)+41611) 3 
9 2187 , 1062882 
1 1 520 2 2(120528C(3)- 16033) 3 

9 ' 2187 531441 


5 


Bose 


A« = 
Ai = 

A 2 = 


1, 14 2 | 768C(3)-311 3 
5 25 ~ l ~ 625 

1, | 2 2 768C(3)+29 3 
5 5 625 
, , 6 r 2 384<(3)+257 3 
fc " r 5 fc 125 


6 


l/Z-cubic 


A« = 
A 2 = 


\ _ 1 38 2 | 2592C(3)-937 3 

Al gt 81 t -|- 21g7 

, 34 2 2592C(3)+1603 3 
fc ' 27 729 


7 


I-tetragonal 


Ai = 
A 2 = 
A 3 = 


1 1 520 2 2(120528C(3)-16033) 3 

9 *~ 2187 531441 

1 250 2 , 5(36936C(3)+41611) 3 
9 2187 1062882 
I 250 2 246888<(3)+165287 3 

t ~ t 243 118098 


8 


II-tetragonal 


Ai = 
A 2 = 

A 3 = 


11 ll c 

11 c ~ 605 

1 58 2 3(123600C(3)+71621) 3 

fc " l ~ 55 fc 166375 
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where e = | corresponds to the physical space dimensionality D — 3. The same 
expansion holds within the plane (u, v). Note, that expression (|T9|) coincides with that 
found for the cubic model with the complex order parameter || . Such a coincidence is 
not occasional because, as was already emphasized, for v — z the model (|l|) goes over 
into to the complex cubic model. 

Unfortunately, RG expansion ([19]) is known to be divergent. Nevertheless, the 
physical information may be extracted from it, provided some resummation method 
is applied. Since the series of N c is alternating, the Borel transformation combined 
with its proper analytical continuation may play a role of such method. To perform 
analytical continuation the Pade approximant only of the type [1/1] may be used within 



given approximation. The Pade-Borel summation of the expansion ([II]) gives 




N c = a + ^rexp )eA — ), (20) 



where a, b, c are the coefficients before e°, e 1 , e 2 in Eq . fll9D , respectively, and Ei(x) is 
the exponential integral. Setting e = ~ in Eq.QSOp, we obtain the value of the critical 
dimensionality 

N c = 1.50. (21) 

This number is close to N c = 1.47 found within the three dimensional RG approach 
[0]. Since N c lies below two, the critical behavior of antiferromagnets (N = 2, N = 3) 
and Nb02 (N = 2) must be governed by the Il-tetragonal fixed point. 

Now let us turn to calculating the critical exponents. To this end, substitute the 
coordinates of fixed points (see Sec. |3|) into the expressions for 7 _1 and r\ (Eqs. (|i"0|) 
and (0)). For the stable fixed point 8 it gives 

r l = l + p^3(l-AT) + p i^(iV-l)(40iV 2 -214iV + 205) 

+ (5jV £ _ 4)5 (l - iV)(12C(3)(5JV - 4)(32iV 3 - 156iV 2 + 159iV - 13) 

- 940iV 4 - 6748iV 3 + 42681 A^ 2 - 67102 A^ + 32558), (22) 



(SA^^" 1 '^- 1 '-^-^ 



(JV-l)(2JV-l) + - 7 — ^ --(iV-l)(190iV 3 -535iV 2 



+ 652A^ - 324) + _ )(1 - A^)(96C(3)(5iV - 4)(32A^ 4 - 128N 3 

+ 212N 2 - 153N + 33) - 10570iV 5 + 22691A^ 4 + 68527iV 3 
- 280399iV 2 + 326888iV - 127676). 
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From Eqs.p 



we find for N = 


2 


I' 1 = 


e 
l -2 


V = 


6 2 

12 + 



7e 2 + e 3 (84C(3) 



24 144 

5e^ e 4 (13-21C(3)) 
36~ + 108 



(23) 



and for N = 3 



_! 6e __ 14e 2 2e 3 (912C(3) + 3905) 

7 ~~ ~TT~T2T + 14641 ' 

10e 2 177e 3 e 4 (50083 - 59328C(3)) 

•n = 1 1 *±m. (24) 

' 121 1331 322102 ' v ' 

where e = | as before. Other critical exponents are found from the well known scaling 
relations. 

We will focus first on qualitative discussion of the results obtained. As was found 
in Ref . , the critical exponents of the Heisenberg and the II-tetragonal fixed points 
coincide within the two-loop approximation. Three-loop analysis yields for the Heisen- 
berg fixed point at iV = 2 

7 -i = 1 _£_7^ + e 3 (28C(3)-ll) 



2 24 48 

e 2 5e 3 e 4 (13-21C(3)) 

V = — + — + — ^- LL - 25 

' 12 36 108 v ; 

Comparing ( pp] ) with (^) we see that the critical exponent 7 of the Heisenberg and 
the II-tetragonal fixed points is different in third order in e, although that difference 
is not too strong. This is one of the results of our investigation. 

It is known that RG series for critical exponents are badly divergent. However 
they contain important physical information which can be extracted provided some 
procedure making them convergent is applied. Unfortunately, it is impossible to use 
simple Pade-Borel summation to process series (|23|)-(|2"5|) because their coefficients have 
irregular signs, in contrast to the critical dimensionality iV c . The most appropriate re- 
summation scheme known for now is a modification of the Borel technique. Principal 
underlying ideas of this method are the analytical continuation of the Borel transform 
beyond its circle of convergence over the cut-plane and a conformal mapping sending 
the cut -plane onto the circle. Such an operation leads to integration of a holomorphic 
function represented by an absolutely convergent series and allows to perform integra- 
tion prior to summation, thus substantiating the perturbation theory approach. The 
algorithm just mentioned incorporates both exactly calculated first several terms and 
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high order asymptotic behavior of perturbation series. For the simple 0(n)-symmetric 
model the coefficients at large order k were shown to look like (— l) k k\a k k b |[I[ [32|. 
It can be expected that in complex models with more than one coupling constants 
asymptotics of RG series will comprise such a factor, at least. Parameters a and b play 
an essential role in the modified Borel method. For a given series 



F(e) = £ f k e k 



transformation 



F (e ~ / e"« - Bit)-, 
Jo V oe / ae 



where S(i) = ^ fc a h r(b+k+i) t k i 1S followed by the conformal mapping 

\/FTI-i 

uj = - — . 

y/t+l + l 

Function B(t) is represented by the series in uo 

B(t(u))=(^- ] j J2A k (X)u k 

where the additional parameter A is introduced to eliminate possible singularity at 
uj = 1. Since the type of that singularity is unknown A is chosen so as to ensure the 
most rapid convergence of the series |28| . 



The main obstruction for application of the method just outlined to the model (|l|) 
as well as to a great deal of other complex anisotropic systems is unknown asymptotic 
parameters a and b. Evaluating them requires enormous efforts. In the case of the 
n-vector model with one coupling constant the parameters a and b have been exactly 
calculated that allowed to obtain accurate numerical estimates for the critical expo- 



nents |24|, |28], [33|, Q. Attempts to find asymptotic parameters for the cubic model also 
were made |35| . They proved to be successful, however, only within the assumption of 
very weak anisotropy. Despite there is no information about asymptotic parameters of 
the model ([!]) available at the moment we chose to resort to the resummation scheme 



of [24, 28 1, in view of the following arguments. Although the asymptotic parameters 



for the isotropic model are explicitly calculated, in Ref. [24| parameter b was varied in 



a neighborhood of the exact value. It is justified by that exact a and b determine large 
order behavior of F(e) while actually one deals with only few terms of perturbation 
series. We believe therefore that, in connection with the model (|I|), similar manipula- 
tions may be valid not only with respect to parameter b but to parameter a as well. 
Variation of a and 6 in a range containing exact asymptotic parameter values of the 
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Table 3: Critical exponents r\ and 7 of the model (1.1) for N = 2 and N = 3 calculated 
within the three-loop approximation 



Type of 
fixed point 


N = 2 


N = 3 


T) 


7 


T) 


7 


Heisenberg 


0.0285 ±0.0002 


1.368 ±0.004 


0.0271 ±0.0002 


1.440 ±0.005 


Bose 


0.0279 ±0.0002 


1.265 ±0.011 


0.0279 ± 0.0002 


1.265 ±0.011 


II-tetragonal 


0.0285 ±0.0002 


1.355 ±0.015 


0.0281 ±0.0002 


1.380 ±0.008 



0(n)-symmetric fixed point and using A as an optimizing parameter result in values of 
the critical exponents displayed in the Table [| Here we suppose that unknown exact 
asymptotic values a and b of the model (|l|) are not much distant from those of the 
O (^-symmetric model. The error of the numerical estimates is established through 
the dispersion of the output due to the variation of a, b, and A. 

As may be seen from the table, the critical exponents of the II-tetragonal fixed 
point appear to be close to those of the Heisenberg point. Unfortunately, we cannot 
compare the critical exponent values of the II-tetragonal fixed point with their two- 
loop analogs and therefore decide how far they shift from the Heisenberg ones with 
higher-loop terms being taken into account. The point is that the estimates of the 
critical exponents were done in Ref. |3| by direct summation of the e-expansion terms 
setting e = 1, that was illegal for the asymptotic series. Under such circumstances, let 
us compare the results obtained with predictions given for the investigated model by 
the RG procedure in three dimensions. Two- and three-loop calculations carried out 
in Refs. [|TJ], [16[ shown that the critical exponents of the II-tetragonal point turned 
out to be close not to those of the 0(2iV)-symmetric model, as in the case of the e- 
expansion method, but to the exponents of the 3D XY (Bose) model. Within the RG 
analysis in 3-D it is a consequence of the closeness of the stable fixed point 8 and the 
Bose point 5 on the three dimensional RG flow diagram. Despite of such a distinction 
in estimates of the critical exponents given by these two RG approaches, one can 
hope that involving higher perturbation orders and using an appropriate resummation 
technique will soften this discrepancy. Not so strong difference between the critical 
exponent values for the II-tetragonal fixed point obtained within the e-expansion and 
the 3D RG methods may serve as a possible confirmation to this conjecture. Indeed, 3D 
RG analysis of the model (|l|) yields the following estimates for the critical exponents 
of the II-tetragonal fixed point: 7 = 1.336, 77 = 0.0261 at iV = 2 and 7 = 1.329, 
7] = 0.0261 at N = 3. Comparing these numbers with their analogs from Table |3|, we 
conclude that the relative deviation does not exceed 4% for 7 and 8% for 77, that is 
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not so bad for the three-loop approximation. An additional stimulus for our hope is 
the beautiful agreement of numerical estimates of the critical exponents for the simple 
0(n)-symmetric model achieved in sufficiently high orders of the perturbation theory 
between the 3D RG [24. 33] and e-expansion [34] approaches. So, for the Ising model 



the magnetic susceptibility exponent was found to be 7 = 1.241 in the frame of 3D RG 
and 7 = 1.239 within the e-expansion method. The relative deviation of these values 
is about 0.1%. 

At last, we would like to emphasize, that although the accuracy of the estimates 
of the critical exponents achieved in the paper cannot be regarded as satisfactory the 
numerical values of the critical exponents for the II-tetragonal fixed point presented 
here are, in our opinion, the most realistic among those so far obtained on the base of 
the e-expansion method. 



5 Conclusion 

The complete RG analysis of a model with three quartic coupling constants and 27V- 
component real order parameter field describing phase transitions in certain cubic and 
tetragonal antiferromagnets as well as the structural phase transition in Nb02 crystal 
has been carried out within the three-loop approximation in D = 4 — 2e dimensions. 
Perturbation expansions for the /3-functions of the record length were obtained using 
dimensional regularization and the minimal subtraction scheme. Coordinates of the 
fixed points and their eigenvalue exponents were calculated for arbitrary N. The 
analysis performed for the eigenvalue exponents has shown that for N > 2 the II- 
tetragonal rather than the Bose fixed point is absolutely stable in the physical space 
within given approximation. The three-loop e-expansion for the critical dimensionality 
of the order parameter N c was found and processed by the Pade-Borel resummation 
technique. The numerical estimate N c = 1.50 obtained confirms the conclusion about 
the stability of the II-tetragonal fixed point. Consequently, the phase transitions in 
the Nb02 crystal and antiferromagnets TbAii2, D1/C2, KiIrCl§, TbD2, and Nd are of 
second order and their critical thermodynamics should be controlled by this point, in 
the frame of given approximation. 

It was observed that the degeneracy of the eigenvalue exponents in the one-loop 
approximation for the II-tetragonal fixed point resulted in certain difficulties in cal- 
culating their e-series. According to the analysis carried out, two-fold degenerate 
eigenvalue exponents should be expanded not in e but in y/e. Although non-integer 
powers of e was shown to drop from the expansions for all N excepting N = 1,2, such 
a degeneracy led to reduction in length of the RG series for eigenvalue exponents and 
therefore to the loss of accuracy expected from given approximation. Indeed, within 
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the three-loop approximation we actually obtain two-loop-like pieces of the series and 
evaluation of the term of order e 3 may require to account the five-loop contributions. To 
understand the structure of the eigenvalue exponent series for the special case N = 2, 
important physically, one has to consider at least four-loop approximation. 

Perturbation expansions for the critical exponents 7 and 77 were calculated up to 
e 3 and e 4 , respectively. For N = 2 the magnetic susceptibility exponents for the II- 
tetragonal and Heisenberg fixed points were found to be different in third order in 
e. For the first time within the e-expansion method the numerical estimates of the 
critical exponents of the model under consideration were given on the base of the 
Borel summation technique modified with a conformal mapping. For the physically 
interesting cases N = 2 and N = 3 the critical exponents of the II-tetragonal fixed 
point turned out to be numerically close to those of the Heisenberg one. On the 
contrary, in the frame of the field-theoretical RG approach in three dimensions the 
critical exponents of the Bose and the unique stable fixed points are close to each 
other. Possibly, these two alternative RG approaches will be in better agreement, 
provided the higher-loop contributions are taken into account. 

The results achieved in our study seem to be self-consistent although there is defi- 
nite discrepancy with the non-perturbative theoretical predictions. We believe that it 
is the effect of insufficiently high approximation employed and the problem of bring- 
ing the results given by the e-expansion method into accordance with those of other 
theoretical approaches and experimental data needs to be solved. 
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A Appendix 

In this appendix we present the eigenvalue exponents of all the fixed points for arbitrary 

N. 

1. Gaussian fixed point 

A u = A„ = \ z = e. 
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0(2N)— symmetric or Heisenberg fixed point 

\u = -e + I1 4 iF (3iV + 7) e 2 -^ i F(48C(3)(iV + 4)(5iV+ll) + 
337V 3 + 2697V 2 + 10727V + 1196V 3 , 



A, = A, = ^L_(JV-2)e+ I ^(57V 2 + 77V + 38)e 2 - 

2pT4yr(48C(3)(iV + 4) (27V 2 + 77V + 23) - 137V 4 - 1997V 3 - 
2(183iV 2 + 98iV-532))e 3 . 

3. Ising fixed point 

\ u = A 2 = |e-fe 2 + ^(2592C(3)-937)e 3 , 
A, = -e + §f e 2 - 7 ig(2592C(3) + 1603)e 3 . 

4. Cubic fixed point 

A i = ~ e + 27^+i) (177V 2 - 27V + 53)e 2 - 

1 ^1296C(3)7V(47V 6 + 47V 4 + 277V 3 + 15iV 2 - 117V - 7) + 



1458^ (Af+l) 3 

32067V 7 - 116837V 6 + 480127V 5 + 345227V 4 - 1118307V 3 + 712057V 2 + 

34527V- 11236) e 3 , 
A 2 = ^e + ^#^(197V 3 - 367V 2 -1657V +106)e 2 + 

4374^(^+i)a (l296C(3)7V(27V 6 + 37V 5 - 227V 4 - 397V 3 - 37V 2 + 207V+ 

7) - 9377V 7 + 78507V 6 - 406747V 5 + 68327V 4 + 1462877V 3 - 996427V 2 - 

271967V + 22472) e 3 , 
A* = ^e + |=P(197V 2 - 1277V +106)6 2 + 

i 3^(l296C(3)7V(47V 4 - 207V 3 + 47V 2 + 217V - 7) - 18747V 5 - 

99977V 4 + 941597V 3 - 1686267V 2 + 1090287V - 22472) e 3 . 

5. Bose fixed point 

\ u = I e -M e 2 + _i_ (768C( 3 ) _3n )e 3 ) 

A x = -I e +| e 2_^(768C(3) + 29) e 3 , 
A 2 = -e + |e 2 - I | g (384C(3) + 257)e 3 . 
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6. VZ cubic fixed point 

X u = A 1 = |e-f e 2 + 11 ^(2592C(3)-937) e 3 , 
A 2 = -e + ff e 2 - yig(2592C(3) + 1603)e 3 . 

7. I-tetragonal fixed point 

A i = ^e + ^(19N 2 -127N+106)e 2 + 

^3^5 (l296C(3)7V(47V 4 - 207V 3 + 47V 2 + 217V - 7)- 

18747V 5 - 99977V 4 + 941597V 3 - 1686267V 2 + 1090287V - 22472V 



A 2 = ^e + ^^(197V 3 -367V 2 - 1657V +106)e 2 + 

4374^(?+i)3 (l296C(3)iV(2iV 6 + 3iV 5 - 227V 4 - 39 iV 3 - 37V 2 + 207V+ 
7) - 9377V 7 + 7850iV 6 - 40674iV 5 + 68327V 4 + 1462877V 3 - 996427V 2 - 
271967V + 22472 V 3 , 



-e+ 2I ? ( -\n (l7N 2 -2N + 53)e 



27N 2 (N+1) [ 

1 ^1296C(3)7V(47V 6 + 47V 4 + 277V 3 + 157V 2 - 117V - 7) + 



1458JV 4 (AT+1)3 

32067V 7 - 116837V 6 + 480127V 5 + 345227V 4 - 1118307V 3 + 712057V 2 + 



34527V -1 1236) e 3 . 
8. II-tetragonal fixed point 

Ai = M e + (sy-i^-i) {4si9n(N - 1)|57V 3 + 67V 2 - 487V + 32 1 - 
3(407V 3 - 2087V 2 + 2537V - 66))e 2 , 



A 2 = M e + (57v-%" ( 2/v-i) {^ign(N - 1)|57V 3 + 67V 2 - 487V + 32| + 

3(407V 3 - 2087V 2 + 2537V - 66))e 2 , 
A 3 = -t + w ^k W zT)(60N 3 -l60N 2 + l8lN-85)e 2 + 

2(5^-4)1(1-2^)3 (48C(3) (27V - 1) 2 (57V - 4)(327V 4 - 1287V 3 + 2127V 2 - 
1537V + 33) + 205607V 7 - 1653287V 6 + 6443927V 5 - 14068647V 4 + 

17567457V 3 - 12243417V 2 + 4337047V - 59052) e 3 . 
Here e = \ corresponds to the physycal space. 
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B Appendix 

As was noted in Sec. [| the II-tetragonal fixed point has an unusual structure of the 
series of the e-expansion for the eigenvalues of the stability matrix. Namely, those 
series are shorter by one order, comparing to their analogs for the other fixed points. 
This is, actually, a consequence of the multiplicity of the roots of the characteristic 
equation in the one-loop approximation, that may cause non-integer powers of e to 
contribute to the expansions. Such a conclusion seems so exotic that deserves thorough 
investigation, to which the present Appendix is devoted. The result of the analysis is 
that for N ^ 1,2 non-integer powers do not appear in the eigenvalue exponents series 
in all orders of the perturbation theory. As to the physically important case N = 2, 
we cannot make such a statement within the three-loop approximation. To answer 
the question whether or not non-integer powers of e will appear in the expansions the 
higher-loop (at least four-loop) approximations need to be considered. 

An eigenvalue A of the stability matrix is a root of its characteristic polynomial. It is 
convenient, rather, to deal with the quantity y — - which is a root of the corresponding 
reduced polynomial denoted hereafter P(y, e). In every order of the perturbation theory 
its coefficients are also polynomials in e, therefore a piece of the series of y(e) determined 
within corresponding approximation coincides with that of some algebraic function. 
Such a function is not analytical in those points on the complex plane where the 
defining polynomial {P{y, e) in our case) has multiple roots. Instead, it has branching 
of an order not greater than the multiplicity of the root. As to the II-tetragonal point, 
at e = (one-loop approximation) the reduced characteristic polynomial has two equal 
roots of the three. It leads to the conclusion that y(e) should be expanded not in e but 
in y/e as a Puiseux series [[36]. That is how half-integer powers of e may occur in the 
series of eigenvalue exponents. Let us show, however, that in the model (H) they are 
absent at least for N ^ 1,2. Consider the reduced characteristic equation 



-y 6 + ay 2 - by + c = (26) 



and assume 



a = clq + a\€ + c^e 2 + 



b = b Q + b l e + b 2 e 2 + ..., (27) 



c = c + cie + c 2 e 2 + . . . , 

V = yo + yie? +y 1 e + yiez + .... 



Here we have omitted higher terms relevant to higher than three loops. The coefficients 
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in ( p7|) are rational functions in TV: 

a ° = (5F34)( 7iV - 8 )' 

ttl = (5JV |_ 4) 3 (270iV 3 - 11297V 2 + 15917V - 736), 

a 2 = 2(5 ~1 4)5 (48C(3)(57V - 4) (1447V 4 - 7207V 3 + 12897V 2 - 947 TV + 230) 

+ 100307V 5 - 1042297V 4 + 4297477V 3 - 8046327V 2 + 6916207V - 222720), 

b o = W ^ W (N-2)(11N-10), 

b x = (5jv ~_ 2 4) 4 (5107V 4 - 31577V 3 + 66157V 2 - 58327V + 1868), 
b 2 = \ (48C(3)(57V - 4) (2727V 5 - 18247V 4 + 44557V 3 - 50957V 2 

+ 27547V - 558) + 258907V 6 - 3384377V 5 + 15470507V 4 - 34371827V 3 

+ 40442037V 2 - 24307527V + 589412), 

c o = ( 5 jv"-4)2 ( N ~ 2 ) 2 ' 

Cl = ^ 4 (TV - 2)(1507V 3 - 8097V 2 + 12297V - 566), 

c 2 = 2(5 ~i 4)6 (48C(3)(7V - 2) (57V - 4) (807V 4 - 4647V 3 + 8657V 2 

- 6417V + 164) + 139507V 6 - 1847457V 5 + 8877057V 4 - 20720607V 3 

+ 25410947V 2 - 15756407V + 389512). 

(28) 

Substituting (|27j ) into ( |26| ) and letting e = we find that y is two-fold degenerate 

taking the value —1 once and 5 2 ^5 4 twice. Comparing factors before equal powers of 

e in (0), we recursively evaluate next terms. The first appearance of a coefficient yi, 

I > 0, occurs at the Z-th step, where it has the multiplier d y P(yo, 0). Factor d y P(yo, 0) 

vanishes due to the multiplicity of yo, hence yi actually cannot be determined from the 

l-th order. So, for y\ from first order in e we have the quadratic equation 

2 oiyg - hyp + d 

y± - — 5 • i^ y j 
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The denominator on the right-hand side is non-zero because it is proportional to 
9yP(yo, 0) and only two of the three roots coincide. Substitution of (p8|) into (|29|) gives 

7/i = 0. The next order (es) does not provide y%, as it might be expected, because at 
this step equation ( po]) vanishes identically. Considering factors before e 2 , we come to 
the quadratic equation for 7/1: 

ViiSyo + oo) + Z/i(2ai7/ - b x ) + a 2 7/ 2 - 6 2 2/o + c 2 = 0. (30) 

The highest coefficient — 37/ + a is proportional to the non-zero quantity 9 2 P(t/ , 0), 
therefore 7/1 is explicitly determined: 

3(7V - 1)(407V 3 - 2087V 2 + 2537V - 66) 
Vl ~ (27V -1) (57V -4) 3 
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4|(iV-l)(iV-2)(iV + 4)(5iV-4)| 

(2iV-l)(5iV-4) 3 ' l J 

Factors before ea in Eq. ( p6|) obey the equation 

By 3 _ = 0, 

where 

B = -6y yi + 2a y 1 + 2a 1 y -b 1 . (32) 

We have come to the crucial point of the consideration. Supposing B ^ we have 
2/3=0. Further, this step can be put into the base of mathematical induction in 

proving disappearance of non-integer powers of e. For m > | we have 
= d y P(y , 0)y m +i + (—62/0 + 2a )yiy m+ i + By m + ..., 

were the terms depending only on y\ with I < m have been suppressed. The first two 

terms turn into zero while B 7^ 0, so the coefficients y m , m > |, can be calculated 

recurrently. Suppose now that m = 2£±i with fc integer. The equation on i/ 2k+i can be 

written in the form 

By2k+i = F(yi,y3,...,y2j^i), 
2 22 2 

where F is some polynomial with zero absolute term in it (because a(e), 6(e), and c(e) 

does not depend upon non-integer powers). Recursively we have y^h+\ = as required. 

2 

Thus, the multiplicity of the roots does not give rise to non-integer powers of e 
unless expression ([32]) vanishes. The important fact is that ([32]) is the derivative of 
the quadratic polynomial ( |3Q|) with respect to y%. It implies that ( Effi) turns into zero 
if and only if y\ is two-fold degenerate. Formula ( pTj) shows that this possibility is 
realized only for N — 1 and N = 2. Concerning these two special cases, one can 
see, that ya is determined from a quadratic equation dyP(yo,0)(y3) 2 = ... arising 
from comparing factors before e 3 in Eq. (p6|) . The right-hand side depends on the 
four-loop contributions. So, if it does not vanish, the expansions would contain non- 
integer powers. Otherwise, ya = and the five-loop approximation gives, in its turn, 
a quadratic equation for y 2 . 

In summary, let us have a look at the structure of the one-loop degenerate eigen- 
value exponent series in general. For every / let d{l) be the order of the expansion of 
the reduced characteristic equation in which y\ is determined. It ranges from Z + 1 to 
21. If d(l) = 21 then y\ is found from a quadratic equation with the non-zero highest 

26 



coefficient dyP(yo,0). We shall say that the solution y(e) splits at the step l s if that 
equation gives two different values of yi s . li y(e) does not split at all, it is convenient to 
assign l s = oo. Let us formulate the resulting theorem as the set of four propositions. 

Theorem 

1. Either the characterictic equation has two equal roots in every order of the per- 
turbation theory or its solution splits at a finite step l s . For every half-integer 
number I from the interval [0, l s ] coefficient y t is determined in the order d{l) = 21. 

2. In the case of finite l s coefficient yi s + m is determined in the order d(l s + m) — 
2l s + m for all m > |. 

3. Coefficients y\ with non-integer numbers I < l s are equal to zero. 

4- Non-integer powers of e contribute to the expansions of the eigenvalue exponents 
if and only if l s is a non-integer number. 

We have demonstrated how the theorem works in the frame of the model under 
consideration, and its full proof will be given elsewhere. 
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